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Abstract 

A unified description of i) classical phase transitions and their remnants in finite systems and 
ii) quantum phase transitions is presented. The ensuing discussion relies on the interplay between, 
on the one hand, the thermodynamic concepts of temperature and specific heat and on the other, 
the quantal ones of coupling strengths in the Hamiltonian. Our considerations are illustrated in 
an exactly solvable model of Plastino and Moszkowski [II Nuovo Cimento 47, 470 (1978)]. 
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I. INTRODUCTION 



In infinite as well as in finite systems a type of phase transition, often referred to as 
a quantum phase transition (qpt), may occur at T=0. Such quantum phase transitions 
differ from classical phase transitions, which can happen only in an infinite systems at Tt^O, 
and generally signal a change in the correlations present in the ground state of the system. 
For an infinite system described by a Hamiltonian, H{X) = Hq -\- XHi, which varies as 
a function of the coupling constant A, the presence of a qpt can easily be understood in 
the following manner Generally the ground state energy is an analytic and monotonic 
function of A. However, if [Hq, Hi] = 0, level crossing may come about and the ground state 
energy is no longer analytic nor monotonic. Although there are other valid mathematical 
reasons that lead to the loss of analyticity[l|, the above simple explanation will suffice for 
our purposes and provides a simple means for defining a qpt in an infinite system. At some 
critical value of the coupling constant, Ac, a new ground state comes to pass. For T > two 
possibilities exist: Ac is an isolated point and the rest the phase diagram is analytic (wrt A) 
or a classical phase transition may occur. In the latter case, for example, for a second order 
phase transition, the free energy is no longer an analytical function of A. As one varies A a 
line of singularities occurs at different temperatures which terminates at T=0 at Ac. This 
provides a simple means of determining Ac, the critical value at which a qpt occurs in an 
infinite system. 

In finite systems a qpt can take place, but strictly speaking classical phase transitions 
can not, since at finite temperatures the partition function and all related quantities are 
analytic. At best only the remnant of a classical phase transition may exist j^. Further- 
more, thermal fluctuations about equilibrium values are large particularly in the region 
where this remnant occurs. For example, studies of their effect on an order parameter have 
concluded that, in atomic nuclei, the super-conducting to normal phase transition is washed 
out0, 5|. However, in spite of these problems a phase diagram has been constructed from 
the remnants in an exactly solvable modelj^ by studying the specific heat, C. 

Clearly information about classical phase transitions or their remnants is contained in 
C. As T — s> 0, however, C 0. In spite of this we will show that it is possible to extract 
information about qpts by studying C in the limit when T — > 0. Only some elementary 
concepts from Information Theory are required. 
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II. FORMALISM 



A. General considerations 

Consider a system whose dynamics is described (at T=0) by the following Hamiltonian 
operator 

H = Ho + \Hi (1) 

where [H(),Hi] = 0. At finite temperatures, the Maximum Entropy Principle of Jaynes6|, |7| 
can be used to determine the appropriate statistical operator, p in the following manner. 
Maximizing the entropy, S{p) = Tr[p\ogp], 



subject to the constraints 



and 



yields 



where 



SpS{p) = (2) 

<H>= Tr[pH] = S (3) 
Tr[p] = 1 (4) 



Z = Tr[e~^^]. (6) 



Generally, in statistical mechanics the coupling constant A is taken to be a constant and 
equation([2]) is used to determine the Lagrange multiplier p. However, in the case of a qpt, 
A is no longer constant and a functional relation between f3 and A may be obtained, using 
equation([3]). 

The specific heat is given by 

c ^ -^\^). (r) 

_ d\ d<H> 

and a necessary and sufficient condition for it to vanish at T = is 

{-^<H>), = (9) 



or equivalently 

Clearly Ac, the critical value of the coupling constant at T=0, can be determined from 
equation(l9]) which clearly indicates that information about the qpt is contained in the specific 
heat. On the other hand C will vanish in this limit if 

I- (") 

for all values of A (see equation ( ITOi) ) . We therefore suggest (and will show) that information 
about a qpt should therefore be contained in the factor {^^^)t=o- Note, however, that 



.d<H>. dE, 



' dx ^^=° dx ■ 

since only the ground state is populated at that temperature. If, indeed as has already been 
pointed out, a qpt occurs at a level crossing then two possibilities exist: 1) a discontinuous 
derivative 

GiX) = (^),=oo(A), (13) 

if does not change sign when passing through Ac, or 2) a null derivative, if does 
change sign when passing through Ac. 

Hence, one has a very nice unified means of identifying both phase transitions and quan- 
tum phase transitions. Furthermore, it is not necessary to begin at finite temperatures to 
find where a qpt takes place . 

For finite systems at finite temperatures (T 7^ 0), C is analytic and structures in ^ 
should be indicative of the remnant of a phase transition. Eq.(l9]) allows one to correctly 
determine the position of the qpt. Alternatively can be used in the manner outlined 
above to determine the position of a qpt (see illustrative graphs in the examples discussed 
below). These two procedures should be equivalent. 



(12) 



III. THE PLASTINO-MOSZKOWSKI MODEL 

This an exactly solvable N-body, SU(2) two-level model Q]. Each level can accommodate 
particles, i.e., is A^— fold degenerate. There are two levels separated by an energy gap E 
occupied by A^ particles. In the model the angular momentum-like operators J^, J^, Jy, Jz, 
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with J {J + 1) = A^(A^ + 2)/4 are used. The Hamiltonian to be here employed reads 

ff = EJ,-e[j2- jf-iV/2], (14) 



a. 



and its eigenstates are usually referred to as Dicke-states [9|]. For convenience we set E = 1 
and 
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J, = (l/2)}_^}_^a+,a,., (15) 

i=l a=l 

with corjespouding expressions for J.. J,. This .s a simple yet nontrivial ease of the Lipkin 
model [10|. For now, we will only discuss the model in the zero-temperature regime. The 
operators appearing in the model Hamiltonian form a commuting set of observables and are 
thus simultaneously diagonalizable. 

The ground state of the unperturbed system (.^ = and at T = 0) is | J, J^) = | y, — y) 
with the eigenenergy Eq = —^N. When the interaction is turned on {C^ ^ 0) and gradually 
becomes stronger, the ground state energy will in general be different from the unperturbed 
system for some critical value of ^ that we will call Ac. This sudden change of the ground 
state energy signifies a quantum phase transition. It should be noted that for a given value of 
A^, there could be more than one critical point. The critical values of the nth transition, i.e., 
Ac at that point, can be found from equation [16] below, provided that Ac > and Ac 7^ 00. 

= iV-(2n-l)- ^^^^ 

A. The iV = 2 problem 

We consider first this simple case, since it can be solved analytically. Here the J = 
1-multiplet for two particles is {Jz = — 1;0;+1}. If we label with the letter % the three 
pertinent Jz— eigenstates one has {Ha = —1; — ^; +1} and {ej} = —1; — ^; 1}, respectively 
and 

Z = e-^ + e^^ + e^ (17) 

with 

If = e^ + ee^^-e-^ (18) 
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Moreover 



Tr[pH] =< E >= Z^^-ef^ - ^e^^ + e-% (19) 

and 

= -[e' + ee'' + e-^] - Z-[-e^ - ^e^^^ + e-^][^]. (20) 
Setting = yields 

= -Z[e^ + ee^^ + e-^] + [e^ + ^e''^ - e"^] [e^ + ^e^^ - e"^], (21) 

i.e., 

[2 cosh (3 + e^^] [2 cosh (3 + ^^e^^] = [2 sinh (3 + ^e^^]"^ (22) 

which is the desired function hnking ^ with (3. Consider now the T = hmit, in which 
(3 oo, cosh/5 e^, sinh/3 — e^. In this hmit fl22|) becomes 

(2e^ + e^«)(2e^ + e'e^«) = (2e^ + ^e^^)' = 4e^^ + ^^e^e"^ + ^'e^/^^, (23) 

entaihng 

(^-1)2 = ^e = l, (24) 

yielding the exact value at which the qpt takes place, as demonstrated in 81]. 
Note, however, one could alternatively start with 

GiO = = ^[-(1 + P)e''Z - i-e^ - ee^^ + e-^)(3e^^.]. (25) 

Requiring 

G(0 = (26) 

one obtains in the limit T ^ 

e^{l-O = (27) 

or 

e = 1! (28) 

which is the exact ,^-value at which the qpt takes place. (Note that at ,^ = 1 G{^) changes 
sign.) In accordance with previous considerations revolving around Eq. (fT3l) . it is clear 
that, at T = 0, the function G above suffers a brutal discontinuity at ^ = = 1, since it is 
"infinite" everywhere except there, where it vanishes. 
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FIG. 1: The lowest three eigenenergies of the = 4 case have been plotted as a function of the 
coupling constant A. There are level crossings at Ac,i = 5 and at Ac, 2 = 1, which is in agreement 
with equation [161 The eigenenergies of the full system are e = ib2,ibl — 3A, — 4A. The solid, 
short-dashed and long-short-dashed line correspond to the eigenenergies e = —2, e = — 1 — 3A and 
e = — 4A, respectively. 

IV. NUMERICAL RESULTS 

Let us now discuss the numerical results for the model given in equation [HI In this 
section we set ^ = A. We will consider the case of four and eight particles, respectively. The 
Hamiltonian is constructed by employing the standard angular momentum matrices in the 
appropriate J-multiplet and is then diagonalized. The resulting 2 J + 1 eigenenergies are in 
general a function of the coupling constant A. This dependence on the coupling constant 
ultimately allows for a level crossing to take place at a critical value of Ac. In figure [1] we 
have shown the subset of eigenenergies that lead to two level crossings (qpt's) in the TV = 4 
particle case. Note that the slope of the ground state energy does not change sign. 

We then construct the canonical partition function Z from the full set of eigenvalues. 
One can now determine the specific heat as given by the two equations [30 
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A. The analogous "specific lieat" 

Once the partition function has been constructed from the eigenvalues of the A^-particle 
Hamiltonian, we are able to form the expectation value of the energy as given by the familiar 
canonical ensemble relation below. 

(29) 

The quantity that will be used to map out the phase diagram of the model, which we 
will call C^A5 is given by the derivative of S, with respect to either (3 or A. In this section 
we will focus our attention on the former case. 



c; = —£ip,\) (30) 



d_ 

w 

A plot of is given in figure [2] for a fixed value oi (3 = 110. The value of (3 was an 
arbitrary choice, in order to demonstrate the following point. At finite temperatures, that 
is when (3 ^ oo, the peaks that are found in figure [2] are a signature of a phase transition 
taking place. They are smoothed out due to finite temperature effects. As the temperature 
is lowered {(3 increases), the peaks move together and become smaller in size. This is shown 
in figure [3l When j3 — > oo, the peaks around each critical point coalesce into a single point, 
namely Ac- This is exactly what one would expect at zero temperature; the phase transition 
takes place where the eigenen-ergies become degenerate. 



B. The analogous "'specific heat"' 

The above investigation of the quantity is one way to characterize the quantum phase 
transitions. It is also possible to investigate the qpt's from another viewpoint. In this section 
we will consider the quantity = -^£{P, A). 

In figure m we have plotted the dependence of C{ on j3 for various values of A. It can be 
seen that if the coupling constant is set in a range corresponding to one particular value of 
the ground state eigenenergy, that at low temperatures tends to the value of the slope of 
the given eigenenergy. For example, when 0<A<|, C^— s>0as/3 becomes large. For that 
range of the coupling constant, the corresponding ground-state eigenvalue is e = —2, which 
of course has a slope of zero. Similarly for ^ < A < 1, ^ —3, which corresponds to the 
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FIG. 2: The quantity has been plotted as a function of the coupling constant A, for the = 4 
particle case, with a fixed value of /3 = 110 (see text for a discussion on this point). There are 
two peaks present in the plot, centered around the two critical points Ac of the system. The peaks 
coalesce into a single point centered at Ac as /3 — > oo (see figure [3]). Everywhere else = 0, in 
agreement with equation [9j 
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FIG. 3: The temperature dependence of Cp in the region of Ac = 1 for the = 4 particle case has 
been plotted. The short-dashed, long-short-dashed and solid peaks correspond to (3 = 70,90, 110 
respectively. One can clearly see that as the temperature is lowered {/3 — > oo), that the peaks 
become smaller in size and narrower in width. In the zero-temperature limit, these peaks would 
coalesce into a single point situated exactly at the location of the quantum phase transition. 
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FIG. 4: as a function of /? for various values of A for the = 4 particle case. The long- 
dashed curves (top two curves) correspond to A = 0.1,0.2; the medium-dashed curves correspond 
to A = 0.5,0.75 (in between the two solid curves); the short-dashed curves (lowest two curves) 
correspond to A = 1.1, 1.2. The solid curves correspond to the the critical values of Ac,n = ^, 1- 
Curves that have same dashing style correspond to the same ground state eigenvalue and in the 
zero-temperature limit tend to the slope of that corresponding eigenvalue. At the critical points, 
picks out the average value of the two slopes from the relevant degenerate eigenvalues. 

slope of the ground state eigenvalue e = —1 — 3A. At the critical values Ac, takes on the 
average value of the slope of the two degenerate eigenenergies involved. In figure [5l we have 
plotted the zero temperature limit of as a function of A. There are two discontinuities in 
the figure, corresponding to the values of Ac where the qpt takes place. The horizontal lines 
in the figure correspond to the slope of the current ground state eigenvalue. 

C. The Plastino-Moszkowski model for N = 8 particles. 

It is also of interest to see if the above methodology works for a larger system. In this case 
the slope of the ground state energy as a function of the coupling constant does not change 
sign. We will briefly summarize the results when the model has N = 8 particles present. 
Using equation [ISl we determine that the critical coupling constants are the following values: 
Ac,n = §5 1, 1- For completeness, the 9 eigenvalues of the system are e = ±4, ±1 — 15A, ±2 — 
12A, ±3 — 7A, — 16A. The quantity is shown in figure [6] and is seen to correctly identify 
where the quantum phase transitions occur. In figure [7] we have plotted in the zero- 
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FIG. 5: as a function of A in the zero-temperature limit for the = 4 particle case. The 
horizontal segments of the plot correspond to the derivatives (with respect to A) of the ground state 



the slope of the ground state eigenvalue e = — 1 — 3A. The discontinuities take place at the critical 



of the two slopes of the relevant degenerate eigenvalues. 

temperature limit as a function of A. As in the = 4-particle case, the discontinuous 
jumps seen in the plot correspond to a quantum phase transition taking place. 

V. CONCLUSIONS 

We have here shown that classical phase transitions and quantum phase transitions can 
be described in a unified fashion. Our treatment has relied heavily on the specific heat and 
is also valid for finite systems where only the remnant of a classical phase transition exists. 
The pertinent considerations were illustrated in an exactly solvable model of Plastino and 
Moszkowski. In particular we have shown that information about qpt's can be obtained 
from the quantity -tt^ and that this equivalent to looking at the zero temperature limit of 



eigenvalue for that particular range of A. For example, 3 < A < 1, CjJ 



3, which corresponds to 



values of the system, viz A^ 



0,1, respectively. At the qpt, the value of is the average value 
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FIG. 6: has been plotted for /3 = 110 (for illustration) in the N = 8 particle case. The peaks are 
centered around the critical coupling constants Ac,n = 7) ■gi 5- Recall that in the zero-temperature 
limit the peaks coalesce into a single point located at the critical points, as has been already shown 
for the = 4 particle case (see figure [3|). Note that we have only plotted the first 3 critical points 
to make the plot clearer; the peak located at Ac = 1 is not shown. 
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FIG. 7: CjJ has been plotted in the zero-temperature limit for the = 8 particle case. The 
horizontal segments correspond to the derivative (with respect to A) of the relevant ground state 
of the system for that particular range of A. The discontinuities take place at the critical values of 
the coupling constant, viz Ac,n = 7> i> 1- the qpt, the value of is the average value of the 
two slopes of the relevant degenerate eigenvalues. 
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the specific heat. 
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